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ABSTRACT 

We measure the weak gravitational leasing shear power spectra and their cross-power 
in two photometric redshift bins from the Canada-France-Hawaii Telescope Tensing 
Survey (CFHTLenS). The measurements are performed directly in multipole space 
in terms of adjustable band powers. For the extraction of the band powers from the 
data we have implemented and extended a quadratic estimator, a maximum likeli¬ 
hood method that allows us to readily take into account irregular survey geometries, 
masks, and varying sampling densities. We find the 68 per cent credible intervals in 
the (Tg-Hin-plane to be marginally consistent with results from Planck for a simple five 
parameter ACDM model. For the projected parameter Ss = cr8(H,„/0.3)°'® we obtain 
a best-fitting value of Ss = 0.768)(Q g 3 g. This constraint is consistent with results from 
other CFHTLenS studies as well as the Dark Energy Survey. Our most conservative 
model, including modifications to the power spectrum due to baryon feedback and 
marginalization over photometric redshift errors, yields an upper limit on the total 
mass of three degenerate massive neutrinos of Emi, < 4.53 eV at 95 per cent credibil¬ 
ity, while a Bayesian model comparison does not favour any model extension beyond 
a simple five parameter ACDM model. Combining the shear likelihood with Planck 
breaks the crs-Hm-degeneracy and yields erg = 0.818 ± 0.013 and = 0.300 ± 0.011 
which is fully consistent with results from Planck alone. 

Key words: cosmology: observations - large-scale structure of Universe - cosmolog¬ 
ical parameters - gravitational leasing: weak. 


1 INTRODUCTION 

The physical nature of the major components of current 
cosmological models is still unknown. Nevertheless, a simple 
six parameter model including dark matter and dark en¬ 
ergy - the A dominated cold dark matter model (ACDM) - 
has been proven very successful in explaining a multitude of 
cosmological observations ranging from the radiation of the 
cosmic microwave background (CMB, e.g. Planck Collabo¬ 
ration et al. 2015a) to supernovae (e.g. Riess et al. 2011) 
and large-scale structure (LSS) probes (e.g. Aubourg et al. 
2014). 

The energy densities of dark matter and dark energy, at 
present, are very well constrained by the aforementioned ob¬ 
servations. The next frontier is pinning down the evolution 
of both dark species, and to observe effects from massive 
neutrinos. One promising probe is the growth of structure 
as inferred from cosmic shear: the (very) weak lensing effect 

* E-mail: fkoehlin@strw.leidenuniv.nl 


due to cosmic large-scale structure bending the light per¬ 
pendicular to the line-of-sight between observer and back¬ 
ground galaxies according to Einstein’s equivalence princi¬ 
ple. The coherent image distortions - the shear - due to the 
gravitational potential of a deflector can only be measured 
statistically, which requires averaging over large numbers of 
sources. Therefore, wide-held surveys covering increasingly 
larger volumes on the sky are required in order to improve 
the precision of the measurements. An analysis of the weak- 
lensing signal as a function of redshift is sensitive to the 
growth of structure, and is thereby indirectly sensitive to 
the expansion rate of the Universe as well as to the cluster¬ 
ing behaviour of various matter species: massive neutrinos, 
dark energy, cold dark matter, etc. 

In order to constrain the dark energy equation-of-state 
and its possible time evolution it is hence crucial to measure 
the cosmic shear signal in different redshift slices (Heymans 
et al. 2013; Benjamin et al. 2013; DES Collaboration et al. 
2015) or directly in 3D (Kitching et al. 2014). 

Massive neutrinos also leave their distinct physical im- 
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prints on the matter power spectrum and hence can be 
probed using weak lensing (e.g. Lesgourgues & Pastor 2006 
and references therein). Theoretically it is straightforward to 
study these features directly in Fourier space, i.e., in terms 
of shear-shear power spectra. Traditionally, lensing analy¬ 
ses employ real-space correlation functions for measuring 
cosmic shear. This introduces further complications in the 
comparison of observations with theory (cf. Section 4.3.2 of 
Planck Collaboration et al. 2015b), because different scales 
are highly correlated. Hence, the signal at very non-linear 
scales requires proper modelling in order to avoid any bias 
in the cosmological parameters. This is generally challeng¬ 
ing due to our limited understanding of the effect of baryons 
on the non-linear matter power spectrum (e.g. Semboloni 
et al. 2011; Semboloni, Hoekstra & Schaye 2013). There¬ 
fore, in this paper we apply a method for extracting the 
data in multipole space and in different redshift bins in 
terms of band powers of the lensing power spectrum. In 
order to achieve this we have implemented and expanded 
the quadratic estimator method originally formulated in the 
context of weak lensing by Hu & White (2001). The Hrst 
applications of this technique to measured shear data were 
presented in Brown et al. (2003) and Heymans et al. (2005) 
using the COMBO-17 and GEMS data sets, respectively. 
More recently, Lin et al. (2012) applied the quadratic esti¬ 
mator technique to shear data measured from the Sloan Dig¬ 
ital Sky Survey (SDSS) Stripe 82. Other recent direct shear 
power spectrum analyses include the Dark Energy Survey 
(DES) (Becker et al. 2015) analysis and the SDSS-FIRST 
cross-power spectrum analysis of Demetroullas & Brown 
(2015). All these studies did not split the power spectrum 
analysis into redshift bins yet and the latter two studies em¬ 
ployed a pseudo-C'(^) power spectrum approach, the other 
major technique for direct power spectrum measurements. 
Alsing et al. (2015) recently presented a hierarchical infer¬ 
ence method that also makes direct use of the shear power 
spectrum. 

In this paper we apply our expanded tomographic ver¬ 
sion of the quadratic estimator to publicly available data 
from the Canada-France-Hawaii Telescope Lensing Survey 
(CFHTLenS, Heymans et al. 2012). CFHTLenS is currently 
the statistically most constraining weak lensing data set and 
covers an area of about 154 square degrees on the sky. The 
data include also photometric redshifts which thus allows us 
to carry out a tomographic analysis. As a further benefit to 
the state-of-the-art data, CFHTLenS has already been used 
before in cosmological analyses (Heymans et al. 2013; Ben¬ 
jamin et al. 2013; Kilbinger et al. 2013; Kitching et al. 2014) 
which enables us to directly cross-check our results with the 
literature. 

The paper is structured as follows: in Section 2 we intro¬ 
duce the weak lensing formalism in terms of power spectra. 
In Section 3 we describe the theory of the quadratic estima¬ 
tor approach and generalize it to include tomography. Sec¬ 
tion 4 provides a brief overview of the CFHTLenS data and 
how to perform shear measurements with it. Before present¬ 
ing the extracted lensing power spectra in Section 6, we test 
and validate the method on mock data in Section 5. From 
the shear power spectra we derive cosmological parameters 
and discuss our results in Section 7. Finally we present our 
conclusions in Section 8. 


2 THEORY 

The deflection of light due to mass is a consequence of Ein¬ 
stein’s principle of equivalence and is termed gravitational 
lensing. One particular case of gravitational lensing is weak 
lensing, the very weak but coherent image distortions of 
background sources due to the gradients of the gravitational 
potential of a deflector in the foreground. 

These image distortions can only be measured in a sta¬ 
tistical sense, given the fact that galaxies are intrinsically 
elliptical, by averaging over large numbers of background 
sources. The resulting correlations in the galaxy shapes can 
be used to study the evolution of all the intervening large- 
scale structure between the sources and the observer, in that 
sense the whole Universe acts as a lens. This particular form 
of weak lensing is called cosmic shear and studied best in 
terms of wide-held surveys covering increasingly more vol¬ 
ume in the sky (cf. Kilbinger 2015 for a recent review). We 
intentionally skip a more basic, mathematical introduction 
of gravitational lensing and weak lensing in particular and 
refer the reader for details on that to the standard literature 
(e.g. Bartelmann & Schneider 2001). 

A wide-held observation of the sky as part of a weak 
lensing survey yields two main observables: the ellipticity 
of galaxies and their (photometric) redshifts. The estimates 
of the ellipticity components ei, 62 at angular positions rii 
can be binned into pixels i = 1 , ..., Npi^ and (photometric) 
redshift bins Zp. The averages of the measured ellipticities 
in each pixel are unbiased estimates of the two components 
of the spin-2 shear held, 7 i(n, Zp) and 72 (n, Zp), which is 
sourced by the convergence held k. In the limit of the hat- 
sky approximation the Fourier decomposition of this held 
can be expressed as 

li{n, Zp) ± i 72 (n, zp) = J 

Zp)±\l5{i,Zp)\ 

where tpi is the angle between the two-dimensional vector t 
and the x-axis. To hrst order for the lensing of density per¬ 
turbations the held /3 vanishes in the absence of any system- 
atics. However, we still want to measure it as a systematic 
test and therefore include it in our notation. The Fourier 
transform of the pixel window function is denoted as W(t.). 
This function can explicitly be written out for square pixels 
of side length Upix in radians as 

W{1) = jo cosipe'j jo sin^J^^ , (2) 

where the zeroth-order spherical Bessel function is dehned 
as jo (a:) = sin(a:)/a:. 

The two-point statistics of the shear held can either be 
expressed in real-space correlation functions or equivalently 
in terms of their Fourier transforms, the shear power spectra. 

Following Hu & White (2001) and expanding the nota¬ 
tion to also include tomographic bins we write out the shear 
correlations between pixels rii and rij in terms of their power 
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spectra as: 

(7iv7ijv) = J [C'®®(^) cos^ 2(pt 

+ sin^ 2(pe 

( 72 iM 72 jv) = J [c®® {£) cos^ 2(pe 

+ sin^ 2ifie 

(7iiM72jv) = J (C™ W - C”™ W) sm4(^f 

+ C™(£)cos4v2r]W""We‘"''”‘-"J^ , (3) 

where we have suppressed the arguments of the shear com¬ 
ponents 7 a (ni, zT) for clarity. 

In the absence of systematic errors and shape noise the 
cosmological signal is contained in the E-modes and their 
power spectrum is equivalent to the convergence power spec¬ 
trum, i.e., C™(£) = and = 0 = C'™(£). Shot 

noise will generate equal power in E- and B-modes. 

The E-mode or convergence power spectra can be pre¬ 
dicted for a given cosmological model: 




r^. 9Ax)9Ax) 

4c" Jo ^ ^HX) 



MxY^) 


( 4 ) 


where x is the radial comoving distance, xa the distance 
to the horizon, a{x) the scale factor, Ps{k-,x) is the three- 
dimensional matter power spectrum, and the angular di¬ 
ameter distance is denoted as /k(x)- Note that we use the 
Limber approximation (Limber 1954) in the equation above 
and the indices p, v run over the tomographic bins. 

The lensing kernels Pn(x) ^ measure for the lensing 
efficiency in each tomographic bin fi and can be written as 

9i^ix) = J dX ix ) 

where n,j(x)dx = P/j( 2 )d 2 is the source redshift distribu¬ 
tion. 


3 QUADRATIC ESTIMATOR 

We summarize here the method originally proposed by Hu & 
White (2001) but make also extensive use of the summary 
provided by Lin et al. (2012). Furthermore, we generalize 
the approach to include tomographic redshift bins. 

We start by assuming that the likelihood of the mea¬ 
sured shear held in terms of band powers B is Gaussian over 
most scales of interest for our analysis, i.e., 

^- (2.)^|C(B)|./^ “'‘[G'^4c(B)]-j]. (6) 

where d denotes the data vector with components 

dfj.ai = Xairii, Zfj) . (7) 

It contains the two components of the measured shear in 
each pixel rii per redshift bin (note that the indices are all 
interchangeable as long as the order is consistent throughout 


the algorithm below). The full covariance matrix C is the 
sum of the cosmological signal and the noise 
The latter includes the contribution from shape and mea¬ 
surement errors. We use the set of equation 3 to build up the 
lensing signal correlation matrix, where we label the shear 
components with indices a, h, pixels with indices i, j, and 
redshift bins with indices fi, v. 

= (7a(ni, 2^)76(nj, z^)) . (8) 


Furthermore, the contribution of shape noise to the signals 
can be encoded in the matrix 


Nu 


dij^ahd [iu 1 


( 9 ) 


where denotes the root-mean-square intrinsic ellipticity 
per ellipticity component for all the galaxies and is the 
effective number of galaxies per pixel i in redshift bin 2 ^." 
Thus we assume that shape noise is neither correlated be¬ 
tween different pixels rii, rij, and shear components 7 ^, 76 , 
nor between different redshift bins 2 ^, 2 ^. This is a well- 
motivated assumption as long as the pixel noise of the de¬ 
tector is uncorrelated. 

We approximate the angular power spectra 
with piece-wise constant band powers B^ep{£) of type 6 G 
(EE, BB, EB) spanning a range of multipoles i within the 
band /3. The index ^ runs only over unique redshift bin cor¬ 
relations. This enables us to write the components of the 
cosmic signal covariance matrix as a linear combination of 
these band powers: 


^sig 


-- YJ, Sc0/3^C(a‘'2) [ 
C,0,/3 


Al 


lep 2(^ + 1) 


nJo{^l')Tab){ij) A~ „ '^i{l')Q(ab)(ij) 


( 10 ) 


The term in brackets in the above equation encodes the ge¬ 
ometry of the shear field including masks and its decompo¬ 
sition in Fourier space. The matrices M,; map the redshift 
bin indices p, v to the unique correlations ( possible between 
those: for riz redshift bins there are only (uz -I-1)/2 unique 
correlations because Zfj, x Zv = Zv x z^i. The explicit expres¬ 
sions for these matrices and the matrices I® and Q® are given 
in Appendix A. 

The best-fitting band powers are determined by 

hnding the cosmic signal C®'® which describes the measured 
shear data the best. For that purpose we use the Newton- 
Raphson method iteratively in order to find the root of 
AC/ABa = 0 (Bond, Jaffe & Knox 1998; Seljak 1998). An 
improved estimate for the band powers Ba is found by eval¬ 
uating the expression 


SBa oc ^ i(P~YAB Tr[{dd^ - C)(C-"DaC-')] , (11) 

B 

where we have introduced now the super-index A for a 
particular index combination {(9P). The matrices Da are 
derivatives of the full covariance matrix with respect to any 
band power combination. We skip here a rigorous definition 
of Da and refer the reader to Appendix A for derivations of 


^ The effective number of galaxies per pixel can be calculated 
using equation 17 multiplied by the area of the pixel Q. 
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these expressions. The elements of the Fisher matrix F can 
be calculated as (Hu & White 2001) 

Fab = I Tr(C-iDAC-iDs) • (12) 

In previous work (cf. Hu & White 2001; Lin et al. 2012), 
the inverse of the Fisher matrix was used as an estimator 
of the covariance between the extracted band powers. We 
refrain from following this approach since the inverse Fisher 
matrix is only an approximation of the true covariance in the 
Gaussian limit. Hence, we decided to estimate the covariance 
of the band powers from mock data instead. We present a 
detailed discussion of this approach in Sec. 5.2. 

For the comparison of the measured band powers to the¬ 
oretical predictions, we have to take into account that each 
measured band power Ba = Bc^ep samples the power spec¬ 
tra with its own window function. This can be computed 
by noting that the expectation value of the band power, 
{Bt^gp), is related to the power spectrum at each wave num¬ 
ber + l)C'c 0 (^)/( 2 ' 7 i') through the band power 

window function W(;gp{£) (Knox 1999; Lin et al. 2012), i.e., 

(13) 

i 

where the sum is calculated for integer multipoles F. The 
elements of the window function matrix can be derived as 
(Lin et al. 2012) 

Wcep{t)= ^ |(F~^){C0/3){x'jA)rxr,A(^) , (14) 

where denotes the inverse of the Fisher matrix (cf. equa¬ 
tion 12). The trace matrix T is defined as 

Ti;ep{£) = Tr(C-iDc9;3C-iD,). (15) 

The derivative denotes the derivative of the full covari¬ 
ance C with respect to the power at a single multipole £. We 
write it out explicitly in Appendix A (cf. equation A19). 

The likelihood-based quadratic estimator automatically 
accounts for any irregularity in the survey geometry or data 
sampling while it still maintains an optimal weighting of the 
data. This is important when dealing with real data because 
it allows for employing sparse sampling techniques and it 
can deal efficiently with (heavily) masked data. The whole 
method and in particular its ability to deal with masks are 
tested extensively in Section 5 before we apply it to data 
from CFHTLenS in Section 6 . 


4 CFHTLENS MEASUREMENTS 

In the following analysis we use the publicly available data^ 
from the leasing analysis of the Canada-France-Hawaii 
Legacy Survey, hereafter referred to as CFHTLenS (Hey- 
mans et al. 2012). The survey consists of four patches (Wl, 
W2, W3, W4) covering a total area of ~ 154deg^. Due to 
stellar haloes or artifacts in the images 19 per cent of the 
area is masked. The leasing data we use in this work are 

^ For the cosmological analysis we employ a range 80 Si ^ Si 2600. 
The lower limit is set by the smallest multipole £ included in the 
analysis and the upper limit must include multipoles £ higher 
than the maximum £ used in the analysis (cf. Section 4). 

^ http://www.cfhtlens.org/astronomers/data-store 


a combination of data processing with THELI (Erben et al. 
2013), shear measurements with lensfit (Miller et al. 2013), 
and photometric redshift measurements with PSF-matched 
photometry (Hildebrandt et al. 2012). A full systematic error 
analysis of the shear measurements in combination with the 
photometric redshifts is presented in Heymans et al. (2012), 
with additional error analyses of the photometric redshift 
measurements presented in Benjamin et al. (2013). One of 
the main results of those extensive systematic tests was the 
rejection of 25 per cent of the CFHTLenS tiles (1 deg^ each) 
for cosmic shear studies. In this work we only use the 75 per 
cent of the tiles which passed the systematic tests as outlined 
in Heymans et al. (2012). Note that this causes considerable 
large scale masking in each patch. 

Photometric redshift measurements have also been ex¬ 
tensively tested (Hildebrandt et al. 2012; Benjamin et al. 
2013) and they were found reliable in the range 0.1 < Zb < 
1.3, where Zb is the peak of the photometric redshift pos¬ 
terior distribution as computed by BPZ (Benitez 2000). In 
our analysis we only use galaxies in this redshift range. 

We compile all tiles associated to a particular 
CFHTLenS patch into a single shear catalogue. Coordinates 
in these catalogues are given in right ascension a and dec¬ 
lination 5 of a spherical coordinate system. We deproject 
these spherical coordinates into flat coordinates via a tan¬ 
gential plane projection. We centre the projection, its tan¬ 
gent point, on the central pointing of each patch. In order 
to measure shears from the ellipticity components ei, 62 as 
measured by lensfit, we first divide the deprojected patch 
into square pixels of side length Upix. We estimate the shear 
components Qa per pixel at position n = (xc,yc) from the 
ellipticity components Ca inside that pixel: 


9a{Xc,yc) 


Wi{ea,i — Ca,i) 


(16) 


where the index i runs over all objects inside the pixel and 
the index a is either 1 or 2 for the two shear and ellipticity 
components, respectively. The weights w are computed dur¬ 
ing the shape measurement with lensfit and they account 
both for the intrinsic shape noise and measurement errors. 
The subscript of the coordinates indicates that the position 
of the average shear is taken to be at the centre of the pixel. 
Note that we assume the galaxies are distributed uniformly 
in the shear pixels. Although this is a simplifying assumption 
we argue that it has only minor effects in the measurement 
considering the general width of the band powers. We define 
distances Vij = I?!; —n^j and angles 93 = arctan {Ay/Ax) be¬ 
tween all pixels i, j which enter eventually in the quadratic 
estimator algorithm (cf. Section 3 and Appendix A). 

In each pixel we apply an average multiplicative cor¬ 
rection (1 -I- m) to the measured shear. This is necessary 
because of noise bias in shear measurements (Melchior & 
Viola 2012; Refregier et al. 2012; Miller et al. 2013). The 
multiplicative correction has been computed from a dedi¬ 
cated suite of image simulation mimicking CFHTLenS data 
(Miller et al. 2013). Moreover, we apply to each measured 
ellipticity an additive correction Ca which is computed from 
all the pass-tiles by requiring that the average ellipticity 
must vanish across the survey as a function of galaxy size 
and signal-to-noise (Heymans et al. 2012). For CFHTLenS 
Cl was found to be zero but for C 2 a correction per object 
has to be applied (Heymans et al. 2012). 
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Table 1. Band power intervals 


Table 2. Effective number densities 


Band No. 

£-range 

range 

Comments 

1 

30-80 

720'-270' 

a), b) 

2 

80-260 

270'-83' 

- 

3 

260-450 

83'-48' 

- 

4 

450-670 

48'-32' 

- 

5 

670-1310 

32'-16l5 

- 

6 

1310-2300 

1615-914 

a) 

7 

2300-5100 

914-412 

a) 


Notes, a) Not used in cosmological analysis, b) No B-mode ex¬ 
tracted. 

The i?-ranges are just an indication and cannot be compared di¬ 
rectly to i?-ranges used in real-space correlation function analyses 
due to the non-trivial functional dependence of these analyses on 
Bessel functions. 


The highest multipole ^pix up to which we want to 
extract band powers employing the qnadratic estimator 
method (cf. Section 3) is on the one hand set by the scales we 
want to investigate because of expected modifications due to 
baryon feedback or massive neutrinos (cf. Section 7.1). On 
the other hand the simplifying assumptions of the algorithm 
such as Gaussianity also limit the maximum ^pix. Hence, we 
only probe into the mildly non-linear regime and consider 
a multipole ^pix ~ 2400 as the maximal physical scale re¬ 
solved. This corresponds to an angular scale of 0.15° = 9' 
and thus sets the pixel size Cpix. We keep parameters fixed 
throughout all CFHTLenS patches such as the side length 
of the shear pixels, (Jpix, measured intrinsic shape noise per 
ellipticity component, = 0.279, and band power inter¬ 
vals. Because the sizes of the CFHTLenS patches are very 
different, the largest distance between shear pixels differs. 
Therefore, we limit our analysis to £fleid ^ 80 (correspond¬ 
ing to an angular separation of pixels of about ~ 4.5°), but 
note that even lower multipoles suffer from more sample 
variance. In summary, the physical scales for our analysis 
are 80 ^ ^ 2300, which corresponds to angular scales 

0.15° SS ^ 4.5°. In total, we choose seven band power in¬ 
tervals enclosing these physical scales as shown in Table 1 
for the E-mode signal extraction. The width of each band 
should at least be two times as wide as ineid in order to min¬ 
imize correlations between the bands (Hu & White 2001). 
The band powers for the B-mode signal extraction are the 
same except that we omit the lowest band power. Note that 
the first band power includes scales below £fieid intentionally 
in order to absorb any DC offsets in the data. The last band 
should Include multipoles above ^pix, because the window 
function of square pixels has a long tail to high multipoles. 
In that sense the enclosing bands are designed to catch noise 
and therefore they are dropped in the cosmological analysis. 
We compute the effective number density of galaxies that 
is used in the lensing analysis and in the creation of mock 
data (cf. Section 5) following the definition of Heymans et al. 
( 2012 ): 


'^eff — ^ 


i(E^ 


(17) 


where Q is the unmasked area used in the analysis and w 
is again the lens fit weight. We show all effective number 
densities per patch and redshift bin in Table 2. 


redshift bin 

Wl 

W2 

W3 

W4 

21 : 0.50 < Zb ^ 0.85 

3.36 

2.80 

3.48 

3.25 

22 : 0.85 < Zb ^ 1.30 

2.86 

2.00 

2.63 

2.22 


Notes. Shown is the effective number density of galaxies (cf. 
equation 17) in arcmin“^ for all four CFHTLenS patches per 
tomographic redshift bin used in this analysis. 


Following the conclusions from Benjamin et al. (2013) 
regarding intrinsic galaxy alignments, which we discuss in 
more detail in Section 6, we define two redshift bins 2 i and 
22 in the ranges 2 i : 0.50 < Zb ^ 0.85 and 22 : 0.85 < Zb ^ 
1.30. These cuts are performed with respect to the peak of 
each galaxy’s photometric redshift distribution Zb ■ For each 
of the two tomographic bins we compute the galaxy redshift 
distribution by summing the posterior photometric redshift 
distribution of all galaxies in the bin, weighted by the lensfit 
weight: 


E, 


Wi 


(18) 


The full galaxy redshift distribution is required in the cal¬ 
culation of the theoretical lensing power spectrum (cf. equa¬ 
tion 5) and it is also needed in the creation of additional 
mock data (cf. Section 5). 


5 METHOD VALIDATION AND 

COVARIANCES 

In order to test and validate the algorithm outlined in Sec¬ 
tion 3 we employ two types of mock data: first we make use 
of the publicly available CFHTLenS Clone’^ (Heymans et al. 
2012) and second we use Gaussian random fields (GRFs). 
This twofold approach is necessary since the multipole scales 
we employ in the cosmological analysis of Section 7 are not 
covered in the CFHTLenS Clone. 

The CFHTLenS Clone is a mock galaxy catalogue that 
consists of 184 independent line-of-sight shear (and conver¬ 
gence) maps with a side length of « 3.58°. These were ex¬ 
tracted via ray-tracing through the TCS simulation suite 
(Harnois-Deraps, Vafaei & Waerbeke 2012) which was pro¬ 
duced with the CUBEP^M N-body code (Harnois-Deraps 
et al. 2013). The CFHTLenS Clone is especially tailored 
to CFHTLenS in terms of the redshift distribution of lens¬ 
ing sources and the noise properties including, for example, 
realistic small scale masks (due to stars etc.). In addition 
to these small scale masks, we randomly mask out three 
non-overlapping tiles of « 1 deg^ each per shear field in or¬ 
der to mimic the effect of the additional ‘bad field’ masks 
also employed in the data. These mask typically 25 per cent 
of the total area of a patch (cf. Section 4 and Heymans 
et al. 2012) and their distribution over a patch does not 
show any systematic preferences. The input cosmology used 
in the creation of the CFHTLenS Clone is WMAP5-like (Ko¬ 
matsu et al. 2009) and summarized in Table 3. Eventually, 
we want to extract scales on the order of several degrees 

^ http://vn90.phas.ubc.ca/jharno/CFHT_Mock_Public/ 
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Table 3. Fiducial cosmology of the CFHTLenS Clone and the 
Gaussian random fields (GRFs) 



Da 

Ob 

h 

tls 

0-8 

Sruiy 

0.279 

0.721 

0.046 

0.701 

0.96 

0.817 

OeV 


Notes. Cosmological parameters used in the creation of the 
CFHTLenS Clone (Heymans et al. 2012) which were also used 
to create the Gaussian random field (GRF) realizations. 

from the data. Kilbinger et al. (2013) showed, however, that 
the power on large scales is significantly underestimated in 
the CFHTLenS Clone. 

In order to also validate the signal extraction on large 
scales, we created 184 Gaussian random field realizations 
(GRFs) of shear fields in two tomographic bins. The fields 
are 20 x 20 deg^ each and generated from convergence power 
spectra that have been computed for the same cosmology 
as the clone, using the measured redshift distributions of 
our two tomographic bins and the modified halof it version 
of Takahashi et al. (2012) for the non-linear contributions 
to the matter power spectrum. Source galaxies are placed 
randomly in the fields with an arbitrary but high enough 
density of 10 arcmin“^ per tomographic bin, and the shears 
are linearly interpolated to these positions. We apply the 
mosaic masks of each GFHTLenS patch to all GRF realiza¬ 
tions in turn, and also apply the patch-specific ’bad field’ 
mask pattern masking about 25 per cent of the total area 
of a GFHTLenS patch. When we compile the actual input 
mock catalogues from the GRF shear fields, we also add 
shape noise by resampling the GRF shear from a Rayleigh- 
distribution with width = 0.279 as measured from the 
data. Furthermore, we randomly sample lensfit weights from 
the corresponding tomographic data catalogues such that 
the effective number densities (cf. equation 17) in the GRF 
mock catalogues match the ones in the data (cf. Table 2). 

The (inverse) Fisher matrices calculated in the 
quadratic estimator algorithm (cf. Section 3) are only an 
approximation of the true (inverse) covariance of the ex¬ 
tracted band powers in the Gaussian limit. In the context of 
a cosmological interpretation of the band powers, however, 
additional non-Gaussian contributions due to the non-linear 
evolution of the underlying matter power spectrum are ex¬ 
pected (cf. Takada & Jain 2009). Hence, we will use our 
mock data also for estimating a more realistic band power 
covariance matrix. 

5.1 Signal extraction validation 

The input cosmology is known for the GRFs and the Clone, 
and we apply a realistic CFHTLenS mask to both sets of 
mock data. We extract the leasing power spectrum using 
the quadratic estimator from the GRFs and the Clone and 
compare it to the input power spectrum. In Fig.l we show 
the residuals between the mean of the extracted band pow¬ 
ers and the predicted band powers for the input cosmology 
for patch W3. The Icr-errors on the mean include the scal¬ 
ing by l/v)/V for N = 184 GRFs for each tomographic bin 
correlation. The binning in multipoles I is the same as the 
one we employ in the final data extraction (cf. Section 4 and 


Table 1). Note that for this test we only extracted E-modes. 
For the calculation of the band power predictions we take 
the convolution with the band window matrices (cf. equa¬ 
tion 14) into account but these are computed for only one 
randomly drawn realization of a GRF. This is due to long 
run-time and we have confirmed for patch W2 that the ran¬ 
domly drawn band window matrix is a fair representation of 
the ensemble (since the noise properties of all GRFs are very 
similar). Fig. 1 demonstrates that the quadratic estimator 
algorithm reproduces the input signal to sufficient accuracy 
and precision, especially given the actual noise level of the 
data (cf. Fig. 4). 

5.2 Band power covariance 

The extracted band powers for each of the 184 shear fields 
from the Clone or 184 GRFs per patch can be used to esti¬ 
mate the run-to-run covariance of the band powers: 

.. "e _ 

Cbw[A,B)^- -^_^^(h;^-Ha)(H^-Hs), (19) 

where is the total number of independent realizations per 
patch, B is the mean of each band power per band over all 
realizations, are the extracted band powers per realiza¬ 
tion, and Ascaie is the scaling factor between each line-of- 
sight clone realization and the actual size of a CFHTLenS 
patch®. The indices A and B denote again the previously 
introduced superindices and run over all bands and redshift 
correlations. 

In order to combine the small-scale covariance esti¬ 
mated from the Clone and the large scale covariance based 
on the GRFs, we stitch both matrices together per patch 
by using the GRF covariance and then replacing all val¬ 
ues associated with a band index for which we want to use 
the Clone covariance. Based on the extensive analysis of the 
Clone and the estimation of covariances from it in Kilbinger 
et al. (2013), we decide to use values from the Clone co- 
variance for multipoles I ^ 670 which corresponds to bands 
5, 6, and 7 (cf. Table 1). Note that bands 7, 6, and 1 are 
not included in any cosmological data analysis though (cf. 
Section 7 and Table 1). 

Due to noise the measured inverse covariance Cb(,i) is 
not an unbiased estimate of the true inverse covariance ma¬ 
trix (Hartlap, Simon & Schneider 2007). In order to 

derive an unbiased estimate of the inverse covariance we 
need to apply a correction derived in Kaufmann (1967) so 
that = aKCB(£')- Assuming a Gaussian distribution of 

the measured band powers B(£), this correction factor is: 


where is the total number of independent mocks, i.e. 
184 in our case, and p is the number of data points used in 

^ Note that Agcaie = 1 in the case of GRFs by construction. For 
the clones we follow (Kilbinger et al. 2013) by matching 90 per 
cent (due to overlapping area between the tiles) of 16 CFHTLenS 
tiles (minus three due to the ’bad field’ masking also employed in 
the clones) into one clone field. The ratio of this number over the 
number of used tiles in one patch (i.e. excluding the ’bad fields’) 
is then l/Agcaie- 
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Figure 1. Residuals between the mean of measured E-mode band powers and predicted band powers for 184 Gaussian random field 
realizations of patch W3. The Icr-error on the mean includes the scaling by l/%/~N for N = 184 measurements. The predicted band 
powers use the known input cosmology (cf. Table 3) and take the convolution with the band window function into account. The residuals 
of each redshift correlation are shown from left to right. 


the analysis. In Section 7 we combine the data of all four 
CFHTLenS patches consisting of four band powers in three 
tomographic power spectra per patch, thus p = 12 for each 
‘patch covariance’. 

We compare the correlation matrix derived from the 
stitched covariance matrix with the correlation matrix based 
on the inverse Fisher matrix which is calculated in the 
quadratic estimator algorithm (cf. equation 12) in Fig. 2 for 
patch W3. The correlation matrices are calculated by nor¬ 
malizing the corresponding covariance matrix with the fac¬ 
tor (Maa with Mas = Csi^i){A,B) or Mas = 

^As- We only include E-mode bands employed later in the 
cosmological analysis in this comparison and find that the 
matrix structure in both approaches is very similar albeit 
with the correlation based on the Fisher estimate being 
smoother, as expected. Finally, we compare both approaches 
in terms of their variance as shown in Fig. 3 for each patch 
individually again only for E-modes used in the cosmologi¬ 
cal analysis. From this comparison we conclude that given 
the noise level in our data the Fisher approach still yields 
compatible error estimates. Nevertheless, we decide to use 
the stitched covariance for our subsequent analysis. This is 
also motivated by the fact that future surveys will yield sig- 
nihcantly improved statistical noise levels and thus require 
a proper covariance estimation beyond the Fisher approach. 


5.3 Computing resources 

We want to comment on the computational requirements for 
our tomographic quadratic estimator approach: the general¬ 
ization of the method to include tomographic redshift bins 
is computationally demanding. The dimension of the covari¬ 
ance matrix defined in equation 8 is set by the size of the 
shear field (times two for the two shear components) and 
the pixel scale. Introducing also two redshift bins increases 
the number of entries in this matrix by a factor of four. 
While this is still efficiently calculated in parallel for smaller 
patches like W2 (« 22.6deg^) and W4 (« 23.3deg^), it be¬ 
comes demanding for patches W3 and W1 (e.g. dim(Cw 2 ) = 



Figure 2. Comparison of correlation matrices for CFHTLenS 
patch W3: the stitched correlation matrix (upper right) is com¬ 
pared to the correlation matrix based on the inverse of the Fisher 
matrix (lower left; cf. equation 12). We show only tomographic E- 
mode bins that enter in the final cosmological likelihood analysis, 
i.e., bins 0 to 3 correspond to 80 ^ £ Si 1310 in the low redshift 
auto-correlation bin, bins 4 to 7 correspond to the same £-range 
in the redshift cross-correlation bin, and bins 8 to 11 correspond 
to the high redshift auto-correlation bin (cf. Tables 1 and 2). 

3076^ versus dim(Cwi) = 9340^) even when exploiting mul¬ 
tiprocessing and optimized libraries such as the Intel® Math 
Kernel Library (MKL®). Nevertheless, the data extraction 
including the calculation of the band window matrices takes 
at most a day on typical cluster machines^. The computa¬ 
tionally most demanding part in our current analysis, how¬ 
ever, is the estimation of the covariance between the band 
powers. This required 184 runs on clones and 184 runs per 
GRF realization per patch. The total runtime for these cal- 


® Version number 11.0.4 

24 cores @2.4 GHz, 256 GB RAM 
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Figure 3. The variances calculated from the stitched covariance 
matrix (solid lines) and the inverse of the Fisher matrix (dashed 
lines) for all four CFHTLenS patches from bottom to top: W1 
(black), W2 (red), W3 (orange), and W4 (cyan). From left to right 
we show the variances in the auto-correlation of the low redshift 
bin, in the cross-correlation between the low and the high redshift 
bin, and in the auto-correlation of the high redshift bin. We limit 
the £-range to the one considered in the cosmological analysis. 


culations was on the order of a month on the same cluster 
configuration for one set of 184 realizations. 

Ongoing and upcoming weak lensing surveys come with 
the advantage of at least an order of magnitude increase in 
survey area compared to CFHTLenS and more regular sur¬ 
vey geometries. Therefore, it will be possible to split these 
surveys into a statistically meaningful number of patches 
still containing scales up to several degrees. This will allow 
for estimating the patch-to-patch covariance directly from 
the data via resampling techniques as an alternative to es¬ 
timating it from mock data alone. However, this approach 
limits the lowest multipole scale to the patch-size and the 
run-to-run covariance will be underestimated at scales close 
to the patch-size. Finally, the rapid advance in terms of num¬ 
ber of cores, clock speed, and internal memory of graphics 
processing units (GPUs) presents a solution to the increase 
in complexity when extending our approach to more redshift 
bins, and/or more band powers, and/or larger contiguous 
patch sizes. The advantage of GPUs lies in their customized 
design to solve linear algebra problems very efficiently and 
massively in parallel which meets exactly the requirements 
of the tomographic quadratic estimator approach. We leave 
an update and porting to GPU programming languages for 
future work. 


6 THE CFHTLENS SHEAR POWER 
SPECTRUM 

For each of the four CFHTLenS patches we extract seven 
E-mode and six B-mode band powers enclosing an interval 
of physically interesting scales of 80 ^ 2300 (cf. Sec¬ 

tion 4 and Table 1). Moreover, we consider two broad mid- 
to high-redshift bins (cf. Table 2) per CFHTLenS patch in 
order to perform a tomographic analysis following Benjamin 
et al. (2013). Doing so, we attempt to decrease the expected 


contamination due to intrinsic galaxy alignments which is 
dominant at low redshifts and high multipoles i. Benjamin 
et al. (2013) concluded that any contamination due to in¬ 
trinsic alignments is at most a few per cent for each red¬ 
shift bin combination. We cross-check this conclusion with 
state-of-the-art intrinsic alignment models constrained by 
recent data from Sifon et al. (2015). For the three intrin¬ 
sic alignment models® employed in there we do not find a 
significant contribution of intrinsic alignments to the cos¬ 
mological signal in any of the redshift bin correlations and 
Gscales employed in our subsequent cosmological analysis. 
Based on these results intrinsic alignments will be ignored 
in the modelling of the signal in our subsequent analysis. 

In Fig. 4 we show the extracted E-mode band powers for 
each tomographic bin. For illustrative purposes we combine 
the band powers extracted from each patch by averaging 
them with inverse variance weights. The errors on the signal 
are estimated from the stitched covariance matrix (cf. Sec¬ 
tion 5.2) whereas the extension of the box in ^-direction is 
just the width of the band. Only bands outside the (grey) 
shaded areas enter in the cosmological analysis (thus we omit 
explicitly the ‘noise catcher’ bands, cf. Section 4 and Ta¬ 
ble 1). Note, however, that for the cosmological likelihood 
analysis we do not use the averaged signals, but instead sum 
the likelihood of each patch as described in Section 7.2. 

We extract E- and B-modes simultaneously. As de¬ 
scribed in Section 2 the cosmological signal is contained in 
the E-modes in the absence of systematic errors. Hence, we 
use the B-mode signal as a systematic cross-check and gener¬ 
ally expect it to be zero within errors. We do not extract the 
EB-modes, which would hint at parity-violation in the data, 
because Kitching et al. (2014) found no evidence for EB- 
modes in the CFHTLenS data. Hence, we decided to only 
include the extraction of B-modes as a non-trivial system¬ 
atic check. We show the extracted B-mode signal per tomo¬ 
graphic bin in Fig. 5. For illustrative purposes we averaged 
the B-mode signal again with inverse variance weights over 
all four CFHTLenS patches. In contrast to the E-modes, 
the lu-errors on the B-modes are derived from the B-mode 
part of the Fisher matrix (cf. equation 12). This is a very 
conservative approach since it will generally underestimate 
the errorbars. The masking in the data might cause leak¬ 
age of E-mode power into B-mode power. In principle, this 
should also be captured by the Fisher matrix but as we ar¬ 
gued in Section 5.2 the Fisher matrix underestimates the 
E-mode error in the intermediate multipole regime due to 
the mildly non-Gaussian intrinsic field. This propagates into 
an underestimated B-mode Fisher-error when compared di¬ 
rectly to a run-to-run B-mode error. However, that does 
not pose a problem as long as we can establish that the B- 
modes are consistent with zero using the underestimated er¬ 
rorbars. We assess the consistency of the B-modes with zero 
via a x^-goodness-of-fit measure and find: Xred(^l) = 1-54, 
xlFW2) = 0.93, X?ed(W"3) = 1.07, and xlAW4) = 0.24 
for 15 degrees of freedom, i.e., including all B-mode bands 
except the last one, which was designed to catch only noise 
due to the long tail of the window function of square pixels 


® These models include intrinsic alignment due to intrinsic ellip- 
ticity correlations, i.e. II, and also intrinsic alignment due to a 
gravitational shear-intrinsic ellipticity correlation, i.e. GI. 
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Figure 4. Measured E-mode band powers in tomographic bins averaged with inverse variance weights over all four CFHTLenS patches 
for illustrative purposes only. From left to right we show the auto-correlation signal of the low redshift bin (blue), the cross-correlation 
signal between the low and the high redshift bin (orange), and the auto-correlation signal of the high redshift bin (red). The low redshift 
bin contains objects with redshifts in the range 0.5 < zi ^ 0.85 and the high redshift bin covers a range 0.85 < 22 ^ 1-3. The Icr-errors 
in the signal are derived from a run-to-run covariance over 184 independent mock data fields (cf. Section 5.2) whereas the extend in 
^-direction is the width of the band. Band powers in the shaded regions (grey) to the left and right of each panel are excluded from 
the cosmological analysis (cf. Fig. 5). The solid line (black) shows the power spectrum for the best fitting, five parameter ACDM model 
derived in the subsequent analysis (cf. Section 7 and Table 4). Note, however, that the band powers are centred at the naive £-bin centre 
and thus the convolution with the band window function is not taken into account in this plot, in contrast to the cosmological analysis. 
We present the E-mode signal for each individual CFHTLenS patch in Appendix B. 



Figure 5. Same as Fig. 4 but for B-mode band powers. Note, however, the different scale (linear) and normalization used here with 
respect to Fig. 4; for reference we also plot the best fitting E-mode power spectrum as solid line (black). We show the measured B-modes 
as (black) dots with Icr-errors derived from the inverse Fisher matrix. Based on these signals we define the shaded regions (grey) to the 
left and right of each panel. E-mode band powers in these regions are excluded from the cosmological analysis (cf. Fig. 4 and see text 
for details). We present the B-mode signal for each individual CFHTLenS patch in Appendix B. 


beyond £pix. However, further tests conducted on the GRF 
mock data show that noise from the last band leaks into the 
second-to-last band depending on the pixel-scale, apix, em¬ 
ployed. This is due to the strong oscillatory behaviour of the 
Fourier-transform of a real-space square pixel (cf. Fig. 2 in 
Hu & White 2001) around ^pix corresponding to Upix. The 
oscillations are amplified if the band is noise-dominated. For 
that reason the B-mode in the second-to-last band appears 
to be more significant than the B-modes in the other bands. 
Removing the second-to-last B-mode band power from the 
X^-goodness-of-fit measure yields the following improved re¬ 
duced x^-values for 12 degrees of freedom: X?ed(^l) = 0.92, 


xlFW2) = 0.80, xlFWS) = 0.61, and x?ed(W"4) = 0.23. 
Hence, we conclude that the B-modes in these remaining 
bands are consistent with zero. Therefore, we only use bands 
2 to 5 in the cosmological analysis of the E-mode signal. 

Following Becker et al. (2015) we define the signal-to- 
noise ratio, S/N, of our band power measurements with re¬ 
spect to the cosmological signal in the mock data from which 
we estimate the covariance: 


S/N = 


^meas ^ g 


di 


-^mock ^mock 

Considering only the band powers used in the cosmological 


( 21 ) 
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analysis (cf. Table 1), we detect a cosmic shear signal in 
W1 at 7.1(7, in W2 at 5.5(7, in W3 at 5.7(7, and in W4 only 
marginally at 2.5(7. Note, however, that the above definition 
of S/N depends on the cosmology employed in the mocks. 
A discrepancy between the mock cosmology and the actual 
cosmology preferred by the data decreases the significance 
in general. 


7 COSMOLOGICAL INFERENCE 


After having extracted the shear power spectrum and hav¬ 
ing derived a more robust estimate of the data covariance, 
we can proceed to the next step: the cosmological inter¬ 
pretation of the tomographic signals, employing a Bayesian 
framework. We estimate cosmological parameters p by sam¬ 
pling the likelihood C(p) with a Monte Carlo Markov Chain 
(MCMC) method. In addition to the parameter estimation 
we also want to compare various model extensions to a base¬ 
line model. 

The Bayesian evidence Z is simply the normalization 
factor of the posterior over the parameters p: 


Z = J d’^pC{p)-K{p) , (22) 

where n denotes the dimensionality of the parameter space 
and 7r(p) is the prior. Since the evidence is the average of 
the likelihood over the prior it automatically implements 
Occam’s razor: a simpler theory with fewer parameters, i.e., 
a more compact parameter space, will have a higher evidence 
than a more complicated one requiring more parameters, 
unless the latter model explains the data significantly better. 
If we wish to decide now between models Mi and Mo, we 
can compare their posterior probabilities given the observed 
data D and define the Bayes factor: 


_ Pr(Mi) 
“ .Zo Pr(Mo) ’ 


(23) 


where Pr(Mi)/Pr(Mo) is the a priori probability ratio 
for the two models, usually set to unity unless there are 
strong (physical) reasons to prefer one model over the other 
a prion. In our subsequent analysis we always assume 
Pr(Mi)/Pr(Mo) = 1. A Bayes factor Kip > 1 implies a 
preference of model Mi over model Mq. Kass & Raftery 
(1995) have proposed a quantitative classification scheme 
for the interpretation of the Bayes factor K (or equivalently 
2 In A). 

Evaluating the usually high-dimensional integral of 
equation 22 is a challenging computational and numeri¬ 
cal task. Here, we employ the nested sampling algorithm 
MULTINEST® (Feroz & Hobson 2008; Feroz, Hobson & 
Bridges 2009; Feroz et al. 2013) via its PYTHON-wrapper 
PYMULTINEST (Buchner et al. 2014) in the framework 
of the cosmological likelihood sampling package MONTE 
PYTHON^® (Audren et al. 2012). 


® Version 3.8 from http: //ccpf orge. cse. rl. ac .uk/gf/project/ 
multinest/ 

Version 2.1.4 from www.moiitepythoii.net 


7.1 Theoretical power spectrum 

In Section 2 we described the calculation of the tomographic 
lensing power spectra (cf. equation 4). These encode the 3D 
matter power spectrum smoothed by tomographic lensing 
kernels (cf. equation 5). For the calculation of the matter 
power spectrum, Ps{k;x), we employ the Boltzmann-code 
CLASS^^ (Bias, Lesgourgues & Tram 2011; Audren & Les- 
gourgues 2011). This already includes the non-linear cor¬ 
rections for which we chose to use the halofit algorithm 
including the recalibrations by Takahashi et al. (2012). Fur¬ 
thermore, CLASS allows us to include (massive) neutrinos 
(Lesgourgues & Tram 2011). The main effect of massive neu¬ 
trinos is a redshift- and scale-dependent reduction of power 
which also propagates into the lensing power spectra 
but is smoothed by the lensing kernels of the corresponding 
tomographic bins (cf. Fig. 6). Over the multipole range of 
interest massive neutrinos lower the lensing power spectrum 
by an almost constant factor. This introduces a degeneracy 
with other cosmological parameters that affect the normal¬ 
ization of the lensing power spectrum. 

We follow Harnois-Deraps et al. (2015) to describe the 
modifications of the power spectrum due to baryon feedback: 


F{k,z) 


Pr‘^{k,z) 
pr^{k,z) ’ 


(24) 


where and denote the power spectra with and 

without baryon feedback, respectively. 

The baryon feedback can be computed from hydrody- 
namical simulations. We use in this work the fitting for¬ 
mula for the baryon feedback derived by Harnois-Deraps 
et al. (2015) using the Overwhelmingly Large Simulations 
(OWLS; Schaye et al. 2010, van Daalen et al. 2011): 


F(k, z) = 1- Abary(Aze^®"'’^ , (25) 

where x = log]^Q(fc/lMpc“'^) and the terms A^, Cz, 
Dz, and Ez are functions of the scale factor a = 1/(1 -|- z) 
which are also dependent on the baryonic feedback model 
(cf. Harnois-Deraps et al. 2015 for the specific functional 
forms and constants). Additionally, we introduce here a gen¬ 
eral free amplitude Abary which we will use as a free parame¬ 
ter to marginalize over while fitting for the cosmological pa¬ 
rameters. In Fig. 6 we show the effect of including baryonic 
feedback on the matter and lensing power spectrum, respec¬ 
tively. In contrast to the effect of massive neutrinos baryon 
feedback causes a significant reduction of power in the lens¬ 
ing power spectrum only at high multipoles. However, this 
is also degenerate with the effect of massive neutrinos on 
these scales. Hence, a proper anchoring of the main cosmo¬ 
logical parameters at low multipoles with high precision is 
paramount if one wants to break degeneracies between all 
these effects. Operating directly in multipole space with re¬ 
spect to both theory and data facilitates the identification 
of distinct features in the power spectra. 


7.2 The shear likelihood 

To compare the measured, tomographic band powers (cf. 
Section 6) to predictions (cf. Section 2), we define 

Version 2.4.3 from www.class-code.net 
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Figure 6. Upper panel: the ratio of modified matter power 
spectra over the dark matter only power spectrum. The dashed 
line (blue) shows the effect of the baryon feedback bias in the AGN 
model from OWLS (Schaye et al. 2010; van Daalen et al. 2011) 
using the implementation by Harnois-Deraps et al. (2015) (cf. 
equation 25). The modifications due to three degenerate massive 
neutrinos with total mass T,mu = 0.18 eV is demonstrated by 
the dash-dotted line (red). The redshift for the power spectrum 
calculation is 2 : = 1.05 corresponding to the median redshift of the 
high redshift bin used in the subsequent analysis (cf. Table 2). 
Lower panel: same as upper panel but for the lensing power 
spectrum of the high redshift bin Z 2 : 0.85 < ^ 1.30 (cf. 

Table 2). 


the shear likelihood as a function of cosmological parameters 

p: 

-2lnC{p) = , (26) 

i a, ^ 

where the index i runs over the four CFHTLenS patches (cf. 
Section 4) and the indices a, P run over the tomographic 
bins. Note that we follow all previous CFHTLenS studies in 
ignoring any covariance between the individual CFHTLenS 
patches. 

The components of the data vector per patch are cal¬ 
culated as 

dL(p) = (fiL-(^?^(p))"'°“), (27) 

where the predicted band powers, {H*(t'))™°‘^®*, depend on 
the cosmological parameters p. They are calculated with 
equations 13, 4, i.e., the band window functions are properly 
taken into account. The inverse of the covariance matrix 
is estimated from a large suite of mock data especially 
tailored to CFHTLenS as described in detail in Section 5.2. 


7.3 Models &: discussion 

In the first part of this cosmological analysis we consider 
the shear likelihood without further combining it with any 
other external cosmological probe. The lensing power spec¬ 
trum is most sensitive to cosmological parameters modifying 
its normalization and slope. Therefore, the normalization of 
the primordial power spectrum, ln(10^°Aa), and the fraction 
of cold dark matter, ficdm, are the primary parameters of 
interest. For an easier comparison of our results with the lit¬ 
erature, we also derive the root-mean-square variance of the 
density field smoothed with the Fourier transform of a top- 
hat filter on a scale R = 8 h“^Mpc in real-space, i.e. us, and 
the total fraction of matter in the Universe, Our base¬ 
line model to which we refer subsequently only as ‘ACDM’ 
includes in addition to these parameters three more free vari¬ 
ables: the Hubble parameter h, the slope of the primordial 
power spectrum Ug, and the fraction of baryonic matter fib. 
The ranges for the flat priors on these parameters are listed 
in Table 4. They follow mostly the ranges employed in the 
CFHTLenS studies by Benjamin et al. (2013) and Heymans 
et al. (2013) in order to assure a fair comparison of our re¬ 
sults with these studies. 

Standard model neutrinos have mass (e.g. Lesgourgues 
& Pastor 2006 and references therein). Hence, we follow 
Planck Collaboration et al. (2015a) in including already two 
massless and one massive neutrino with the (hxed) mini¬ 
mal mass of Em„ = 0.06 eV (assuming a normal mass hi¬ 
erarchy with one dominant mass eigenstate) in our baseline 
ACDM model. Moreover, we always assume a flat cosmolog¬ 
ical model. 

The first extension of the baseline model is to introduce 
a free total mass for three degenerate massive neutri¬ 

nos. We refer to this model as ‘ACDM-|-i/’. Since we expect 
the effect of massive neutrinos to be degenerate with the ef¬ 
fect of baryonic feedback, especially at high multipoles (cf. 
Section 2 and Fig. 6) we investigate this effect in the model 
‘ACDMa’: here, we additionally include the fiducial baryon 
feedback model of equation 25 with Abary = 1 for the AGN 
model taken from the OWLS project (Schaye et al. 2010; 
van Daalen et al. 2011). The degeneracy between baryonic 
feedback and massive neutrinos is investigated in the model 
ACDMa-|-i/, where Emi, is free to vary but which includes 
the fixed hducial baryon feedback model. We relax the as¬ 
sumption of a fixed baryon feedback model in the model 
‘ACDM-|-Abary’ by allowing the amplitude of the feedback 
Abary to Vary (cf. equation 25). Combining the assumption 
of a free amplitude in the baryon feedback model and a free 
total mass of three degenerate massive neutrinos, Em,^ in 
the model ‘ACDM + u + Abary’ yields a maximally degener¬ 
ate model in baryonic feedback and neutrinos. In total this 
model consists of seven free parameters. 

Moreover, we want to test the effect of a photometric 
redshift bias which causes a coherent shift of the photometric 
redshift distributions per tomographic bin (cf. equation 5) 
by AZfi. Hildebrandt et al. (2012) showed that the bias on 
photometric redshifts in CFHTLenS is Az < 0.02 (cf. their 
Fig. 8). However, this estimate does not account for out¬ 
liers which can increase the photometric redshift bias sig- 
nihcantly. Therefore, we make a more conservative assump¬ 
tion and treat the photometric redshift biases within 
a flat prior range of —0.05 ^ AZfj, ^ 0.05 as nuisance pa- 
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rameters to marginalize over. In the most complex model 
ACDM + zz + Abary + , which we abbreviate subsequently 

to ‘ACDM + air, we include a free amplitude for the baryon 
feedback model, massive neutrinos and treat the photomet¬ 
ric redshift biases as nuisance parameters. 

All models, their prior ranges and the parameter esti¬ 
mates derived from the likelihood sampling are summarized 
in Table 4, where we always quote the weighted median value 
for each varied parameter. The errors denote the 68 per cent 
credible interval of the posterior distribution after marginal¬ 
ization over all other free parameters. 

We compare the 68 and 95 per cent credible inter¬ 
vals for the baseline ACDM model and the most complex 
ACDM -I- all model in Fig. 7. Both models are marginally 
consistent with the 68 per cent credible interval from Planck 
Collaboration et al. (2015a) (TT+lowP) and the most com¬ 
plex ACDM -I- all model is fully consistent with Planck at 95 
per cent credibility. This model is very conservative because 
it also accounts for a possible photometric redshift bias per 
tomographic bin, and thus is expected to yield the largest er- 
rorbars. For this model we show marginalized ID posteriors 
for every free parameter (cf. Table 4) and marginalized 2D 
contours for every parameter combination in Fig. 8. From 
this figure but also from Table 4 it is apparent that our 
parameter constraints are weaker than those derived from 
Planck. The shear data are also unable to constrain the slope 
of the primordial power spectrum, rig, especially once the 
models also include massive neutrinos, since both parame¬ 
ters influence the slope of the lensing power spectrum in a 
similar way. Hence, the estimate on rig is following the flat 
prior distribution. From our most conservative model ex¬ 
tension, ACDM+all, we derive an upper limit on the total 
mass of three degenerate massive neutrinos at 95 per cent 
credibility of Emi, < 4.53 eV. In contrast, Planck Collabora¬ 
tion et al. (2015a) (TT+lowP) derive an upper limit (95 per 
cent) on the total mass of three degenerate massive neutri¬ 
nos of Em,y < 0.72 eV. Combining the primary CMB data 
with secondary data and/or other external probes lowers the 
upper limit to < 0.17 eV. 

In the (jg-Om-plane we can directly compare to the re¬ 
sults from the CFHTLenS analysis by Heymans et al. (2013). 
They employed a 6-bin tomographic real-space correlation 
approach and in Fig. 9 we show the 68 per cent credible 
intervals for their conservative model including a marginal¬ 
ization over intrinsic alignments. The 68 per cent credible 
intervals of our baseline ACDM model is consistent with the 
one derived by Heymans et al. (2013). However, the contours 
of our model are generally broader because we use only two 
tomographic bins. 

The shear power spectrum is most sensitive to the pa¬ 
rameters Dcdm and ln(10^°Ag) or equivalently to Dm and 
erg. However, as can be seen in, for example, Fig. 7 the re¬ 
lation between Dm and ng is degenerate and what lensing 
can actually constrain best is the combination of both pa¬ 
rameters in the projected quantity (jg(Dm/0.3)“. The value 
of a depends on the scales probed and is connected to the 
width of the likelihood contour. We derive it from fitting the 
function In (Tg(Dm) = —a In Dm+const. to the likelihood sur¬ 
face in the crg-Dm-plane. Since we find it to be consistent 
with « 0.5 in all our models, we follow DES Collaboration 
et al. (2015) in defining the quantity Sg = cr8(Dm/0.3)°'®. 
We present values for this parameter combination obtained 



Figure 7. Shown are 68 and 95 per cent credible intervals (blue, 
inner and outer contour, respectively) for our baseline ACDM 
model where Dmi n^g, h, n^, and Qi, are free to vary. Additionally 
shown are the 68 and 95 per cent credible intervals (red, respec¬ 
tively) for our most complex model ACDM + all where also the 
total mass of neutrinos Sm^, the amplitude for the Baryon feed¬ 
back model Abary, and a systematic photometric redshift bias 
per tomographic bin are free to vary. We marginalize over 

all other free parameters. Finally, we plot the 68 and 95 per cent 
credible intervals derived from Planck Collaboration et al. (2015a) 
(TT+lowP). 



Figure 9. We show the 68 per cent credible interval (blue) for 
our baseline ACDM model. Additionally shown is the 68 per cent 
credible interval for the 6-bin tomographic real space analysis 
from Heymans et al. (2013) (cf. also their figure 4) where intrinsic 
alignments are marginalized over (red). Finally, we plot the 68 per 
cent credible interval from Planck Collaboration et al. (2015a) 
(TT+lowP). 

from the above shear-only likelihood sampling in Table 5. 
We compare the values of Sg for all our models in Fig. 10, 
where we also show the constraint on that parameter combi¬ 
nation by Planck Collaboration et al. (2015a) (TT+lowP). 
For this combination all our tested models are consistent 
with each other. However, all models are in mild tension 
with the constraint on Sg derived from Planck (TT+lowP). 

Moreover, we present in Fig. 10 the constraints on Sg of 
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Table 4. Cosmological parameters from shear likelihood only 


Model 

Ocdm 

ln{10'"Aa) 

Dm 

O’s 

Db 

Tig 

h 

(eV) 

Abary 

A2I 

A22 

Prior ranges 

lo.,i.] 

(0.,10.] 

derived 

derived 

[0.,0.1] 

[0.7.1.3] 

[0.4,1.] 

[0.06,6.] 

[0.. 10.] 

[-0.05,0.05] 

[-0.05,0.05] 

ACDM 

0.211^5 

o Eq+1-49 
'5-0o_i 52 

0.264:2 

0.844:2 

0.054:2 

1 01 +0-29 
1-1^1-0.23 

0 62+^ *^® 
U.DZ_o 22 

- 0.06 

- 

- 

- 

ACDMa 


3.504;« 

0.254:2 

0.854:4 

0.054:2 

1 00+°'^® 
l.UU_9 22 


- 0.06 

= 1. 

- 

- 

ACDM + :/ 

n 21 

'^■'^-‘--0.13 

3.654-2 

0.304-2 

0.754:2 

0 04+'^''^^ 

1 O'y+O-^^ 

i.UO_Q 28 

0 70+0 1® 

i.37i;:?; 

- 

- 

- 

ACDMa + u 

0 21 

3.694:2 

0-294:2 

0.764:2 

0.044:2 

l.UO_Q 27 

0 71 +0-22 

i.34i;:?; 

= 1. 

- 

_ 

ACDM + Abary 


3.624-4 

0.264:2 

0.854:2 

0.054:2 

1 00+°'^® 
l.UU_Q 24 

0.6012.2 

= 0.06 

2-9022?) 

- 

- 

ACDM + 17 -|- Abary 

fv 99+0.08 

'^■^^-0.13 

3.6914) 

0.304:2 

0 76+0 15 


1.06122 

0-694-4 

i.29i;:?) 

2.5144 

- 

- 

ACDM + As;, 

0 94+0-1*3 

q 2fi+l-28 
o.4;o_i 32 

n 2Q+010 

n sn+O-21 

U.8U_o.22 

u.ua_Q 03 

0.98122 

0.624-2; 

= 0.06 

- 

0.0312?; 

—0 02+*^-^^ 
D.U^_0 03 

ACDM + all 

0 24+0-09 

q r'T+1.34 

n q2+010 

U.OZ_Q 13 

0 74+^-^^ 

0 04+0-02 

1-04222 

0.6712;? 

1 q9 + U.Ot) 

i.OZ_l 26 

2.4922;? 

0.031?:?; 

—0 02+'^-*^^ 
U.U.i_0.o3 


Notes. We quote weighted median values for each varied parameter and derive Itr-errors using the 68 per cent credible interval of the 

marginalized posterior distribution. 


Table 5. Constraints on Ss and (Tg(Om/0.3)“ 


Model 

58 = (T8(n„,/0.3)0-5 

Mean error on Ss 

(7g(nm/0.3)“ 

a 

Shear likelihood only 





ACDM 

0.7684°;°« 

0.042 

U./DZ_Q Q 3 g 

0.538 

ACDMa 

0-7704°o:°^^ 

0.043 

n 765+0 044 

U. /DO_q Q 3 g 

0.533 

ACDM -I- u 

n 7'Q7'+0.057 
<3'-0.054 

0.056 

n 707 + 0.057 
^+'^'-0.055 

0.479 

ACDMa -|- u 

0 741+0-0^5 

IJ- '^J-_0.047 

0.051 

0 741+0-056 

O'-' ^1-0.046 

0.465 

ACDM -|- A)3a^j.y 

n 777+0.048 
‘ -0.040 

0.044 

n 770+0.046 
ll-' ''^-0.040 

0.531 

ACDM ly Abary 

n 740 + 0 .055 
U.74S_q Q 4 g 

0.052 

n 748+0-054 

U.745_Q Q 5 Q 

0.479 

ACDM -1- Azf, 

n 77’! +0.050 
' * -^-0.039 

0.045 

0 767+0-045 

0.555 

ACDM -1- all 

n 7 c;c:+0-059 

U. (OO_0 059 

0.059 

n 755+0-059 

U. / 00 _Q Q 5 g 

0.491 


Notes. We quote median values for the constraints on Ss = a'g(flni/0.3)® ® and (Tg(flm/0.3)“. The errors denote the 68 per cent credible 
interval derived from the marginalized posterior distribution. 


other lensing studies. In particular, we compare to the recent 
constraint from DES Collaboration et al. (2015) (“Fiducial 
DES SV cosmic shear”). This study employed a real-space 
correlation function approach in three tomographic bins. We 
find our constraints to be consistent with theirs which is 
mainly due to the large errorbars of the measurement on 5g 
from the Dark Energy Survey (DES). In addition to their 
own results DES Collaboration et al. (2015) also resampled 
the likelihoods of the CFHTLenS studies from Kilbinger 
et al. (2013) and Heymans et al. (2013) and derived con¬ 
straints on Ss- We show these constraints also in Fig. 10. Kil¬ 
binger et al. (2013) employed a non-tomographic real-space 
correlation function approach and their constraint in Fig. 10 
employs “all scales” out to large angular scales U ~ 350'. The 
constraint from Heymans et al. (2013) in Fig. 10 uses only 
the “original conservative scales”. Our results are consistent 
with both these studies, as was already the case for Heymans 
et al. (2013) in the full crg-Dni-plane (cf. Fig. 9). 

For the comparison of our results to other CFHTLenS 
studies and the originally published constraints from Kil¬ 
binger et al. (2013) and Heymans et al. (2013) we have to 
resort to the parameter combination (Tg(Dni/0.3)“. The ex¬ 
ponent a is in general quite similar between the quoted lens¬ 
ing studies but not the same, which the reader should bear 
in mind when looking at Fig. 11. For completeness, we show 
again the constraints on that parameter combination from 
Heymans et al. (2013) and Kilbinger et al. (2013). Kitching 
et al. (2014) employed a 3D lensing approach which allows 


Planck 2015 
DES 2015 
Kilbinger 2013 
Heymans 2013 
ACDM+all 
ACDM + v + A 




ACDM+A, 
ACDMa+n 
ACDMa 
ACDM + Azj, 
ACDM+n 
ACDM 


_L 


h 


J_ 


H - 


X 


0.65 0.70 


0.75 


0.80 0.85 


0.90 


,/0.3) 


0.5 


Figure 10. Shown are Icr-constraints on the parameter combi¬ 
nation = (T8(r2m/0.3)^'^ for all of our tested models (cf. Ta¬ 
bles 4 and 5). We compare them to constraints from other lensing 
analyses and to the constraint from Planck Collaboration et al. 
(2015a) (TT+lowP). Note that for Heymans et al. (2013) and 
Kilbinger et al. (2013) we quote the values derived in DES Col¬ 
laboration et al. (2015) for the “original conservative scales” and 
for “all scales”, respectively. ‘DES 2015’ refers to the fiducial result 
from DES Collaboration et al. (2015) (“Fiducial DES SV cosmic 
shear”). 
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Figure 8. Shown are all parameter constraints from sampling the likelihood of model ACDM + all. The dashed lines in the marginalized 
ID posteriors denote the weighted median and the 68 per cent credible interval (cf. Table 4). The contours in each 2D likelihood contour 
subplot are 68 and 95 per cent credible intervals smoothed with a Gaussian. 


for control over the fc-scales included in the analysis. How¬ 
ever, their constraint on (T8(flm/0.3)“ for which we quote 
the value including large scales, i.e., k ^ 5/iMpc“^, yields 
by far the largest errorbars due to which their constraint is 
consistent with all other CFHTLenS studies and also consis¬ 
tent with the Planck constraint. The analysis by Benjamin 
et al. (2013) is the most similar to the one presented here: al¬ 
though their analysis employed a real-space correlation func¬ 
tion approach and did not include scales as large as the ones 
used here, the two redshift bins in their tomographic analy¬ 


sis are exactly the same ones employed in this analysis. The 
constraints are also consistent with each other and especially 
our ACDM model also yields comparable errorbars. 

In summary, all models are consistent with each other 
mainly due to increasing errorbars for increasingly more 
free parameters. For the comparison of our analysis to 
other cosmic shear studies we derived a constraint on the 
projected parameters Sg or (j8(Hni/0.3)“. In general, we 
find consistency in these projected parameters with all 
other CFHTLenS studies and DES. Employing the Bayesian 
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Table 6. Evidences from shear likelihood only 


Model 

InZ 

2 In A {K = ZilZACDM) 

interpretation 

ACDM 

-40.96 ±0.06 

0 

- 

ACDMa 

-41.07 ±0.06 

-0.22 

support for ACDM 

ACDM -1- n 

-41.63 ±0.07 

-1.34 

support for ACDM 

ACDMa -h 1 / 

-41.83 ±0.07 

-1.74 

support for ACDM 

ACDM -h 

-41.66 ±0.06 

-1.40 

support for ACDM 

ACDM -|- z/ + 

-42.48 ±0.07 

-3.04 

support for ACDM 

ACDM -1- Azi, 

-40.75 ±0.07 

0.42 

preference over ACDM “not worth more than a bare mention” 

ACDM -1- all 

-42.19 ±0.07 

-2.46 

support for ACDM 


Notes. For each model we calculate the global log-evidence, In^, and compare all evidences in terms of the Bayes factor K (or equivalently 
21ni^) to the baseline ACDM model. The interpretation of the Bayes factor is following the scheme proposed by Kass & Raftery (1995). 


Planck 2015 

-1-r 

^ — 1-1-^-1-1- 

hH 

DBS 2015 

t— i-H 

Kitching 2014b 


1-*- 

Kilbinger 2013b 

- h- * —1 

1 

Heymans 2013 

HH 

1 

1 

Benjamin 2013 

- H*H 


ACDM-kail 

- \- 9 -\ 

1 

ACDM 

_1_i_ 

J_ 1 _ 1 _ 1 _ 1 _ 1 _ 


0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 




Figure 11. Shown are Icr-constraints on the parameter combina¬ 
tion (T8(f2in/0.3)“ for our ACDM and ACDM-|-all models (cf. Ta¬ 
ble 4). We compare them to constraints from other lensing analy¬ 
ses and to the constraint from Planck Collaboration et al. (2015a) 
(TT+lowP). Note that for Heymans et al. (2013) we quote the 
value derived by marginalising over intrinsic alignments. For Kil- 
binger et al. (2013) and Kitching et al. (2014) we cite values 
including the largest scales in their analyses. 


model comparison framework, we can decide which of the 
tested models describes the shear data best: in Table 6 
we present the natural logarithms of the evidence for each 
model. Comparing these models in terms of their Bayes 
factor K with respect to the simplest models ACDM or 
ACDMa, we find no evidence for any of the tested extensions 
except for a very weak preference of the model ACDM-I-A^^ 
over our baseline model which is according to the interpreta¬ 
tion scheme of Kass & Raftery (1995) “not worth more than 
a bare mention”. 

Hence, we conclude that the extracted band powers of 
the tomographic shear power spectra measured over a range 
80 ^ 1310 are described sufficiently well within their 

errors by a standard, five parameter ACDM model. 

Finally, we combine our shear likelihood with the most 
recent data and likelihood release'^^ from Planck Collabo- 


PLC-2.0 from http://pla.esac.esa.int/pla/ 


ration et al. (2015a) in order to break the degeneracy be¬ 
tween the parameters Dm and erg. In particular we employ 
the Planck primary CMB temperature data (TT) from high 
multipoles I in combination with the Planck low multipole 
polarization data (lowP). Due to long run-time we chose to 
use the PLIK HIGHL-LITE likelihood code which requires 
only to marginalize over one nuisance parameter, Apianck. 
The Bayesian model comparison showed no evidence for 
any model extension beyond a baseline ACDM model for 
describing the shear likelihood. This implies that we would 
essentially reproduce Planck-orAj results if we were to add 
parameters for which there is no evidence. Hence, we con¬ 
sider only six cosmological parameters and one nuisance 
parameter for the combined ‘Plancfc+Shear’ model: Dcdm, 
ln(10^°Aa), h, Db, ris, Treio, and Apianck- Again we assume 
one dominant neutrino mass eigenstate in the normal hier¬ 
archy with Em,y = 0.06 eV and a flat cosmology. In com¬ 
parison to our shear-only likelihood analysis we chose to 
use narrower prior ranges for most parameters (cf. Table 7). 
Due to the reduced set of nuisance parameters with respect 
to the original Planck analysis, we also resample the Planck 
likelihood for the seven parameter baseline model so that 
comparisons of likelihood contours are fair. 

Prior ranges and parameter constraints for the 
resampled Planck likelihood and the combination of 
Planck+Sheai are presented in Table 7. Fig. 12 demon¬ 
strates that combining the shear likelihood with the Planck 
likelihood yields Improved constraints on as and Dm and 
breaks the degeneracy between the two parameters. The 
68 and 95 per cent credible intervals are largely overlap¬ 
ping and show marginal consistency between the two data 
sets as already observed above. We find the constraints 
as = 0.818 ± 0.013 and Dm = 0.300 ± 0.011 which are con¬ 
sistent with the constraints from the resampled Planck-only 
likelihood (cf. Table 7). 


8 CONCLUSIONS 

In this work we generalized the original quadratic estima¬ 
tor approach by Hu & White (2001) to include tomographic 
bins. We validated the method and its extension to tomo¬ 
graphic bins by applying it to mock data tailored to the 
survey specifications of CFHTLenS. In particular we made 
use of the official CFHTLenS Clone but produced also our 
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Table 7. Cosmological parameters from a combined analysis of the shear and Planck likelihoods 


Model 

^edm 

ln(10i“As) 

Hni 

0-8 

Ob 

Hs 

h 

Tpeio 

^Planck 

Prior ranges 

[0.1, 0.4] 

[2.. 4.] 

derived 

derived 

[0.,0.1] 

[0.8,1.2] 

]0.5,0.8] 

]0.04, 0.12] 

]90.,110.] 

Planck (TT+lowP) 
Planck-\-She&T 

u.z,uo_Q Q]^3 

U.ZOi_Q QJ^Q 

o.uyo_Q Q34 
o /-i'7’7+0.037 

o-'ai '-0,035 

n "^10+0.013 

O.OTO_q 014 

U.OUU_o oil 

0.830l);;S“ 

0.818l«;S)^ 

o o 

b b 

00 CO 

1 +1 + 

o oo o 

b bb b 

o oo o 

U.yDD_o 006 

0 97i+0'006 

u.y (T_o 006 

U.D/'l_o,oio 

n fiR/fTO.oos 

0.d84_o,oo9 

0 079+0-01* 
u.u'y_o,oi9 

U.U'«i_o,oi8 

100.041“;^,^ 

100.021“:^^ 


Notes. We quote weighted median values for each varied parameter and derive Itr-errors using the 68 per cent credible interval of the 
marginalized posterior distribution. For the model Planck (TT+lowP) we resampled a simplified version of the original likelihood that 

includes only one additional nuisance parameter, Apianck- 


0.95 

0.90 

0.85 

0.80 

0.75 

0.26 0.28 0.30 0.32 0.34 0.36 0.38 0.40 


- Planck (resampled) 

- P/ancfc+Shear 



Figure 12. We show 68 and 95 per cent credible intervals (red, 
inner and outer contour, respectively) derived from sampling only 
the Planck likelihood (TT+lowP) with the simplified model con¬ 
sisting of six cosmological parameters and only one nuisance pa¬ 
rameter (cf. Table 7). We combine the Planck likelihood then with 
the shear likelihood and sample from the combined likelihood for 
the same simplified model and derive 68 and 95 per cent credible 
intervals (blue, inner and outer contour, respectively). 


own sets of Gaussian random field realizations in order to 
test the performance for the larger scales used in our analy¬ 
sis. We also used the 184 independent shear catalogues from 
the CFHTLenS Clone and our GRFs to derive a run-to-run 
covariance. We applied the method to public shear data from 
CFHTLenS in two tomographic bins to extract band powers 
of the lensing power spectrum. 

We use the extracted band powers and the run-to-run 
covariance estimated from our suite of mock data in or¬ 
der to sample the shear likelihood. The sampling is per¬ 
formed in a Bayesian framework. We derive constraints on 
cosmological parameters as well as the Bayesian evidence 
for each model. In addition to the five baseline cosmologi¬ 
cal parameters, our most conservative model extension in¬ 
cludes a free total mass of three degenerate massive neutri¬ 
nos, a free amplitude for the baryon feedback model of the 
matter power spectrum and photometric redshift biases to 
marginalize over. For this model we derive an upper limit 
on the total mass of three degenerate massive neutrinos of 
Emi, < 4.53 eV at 95 per cent credibility. Based on the anal¬ 
ysis of the shear likelihood we find no evidence for any of the 
tested model extensions though: a standard, five parameter 
ACDM model is sufficient to describe the lensing power spec¬ 


trum band powers measured over a range of 80 1310 

in two tomographic bins. The main parameters constrained 
by the lensing power spectra are as and flm and we find 
the 68 percent credible intervals in this parameter plane 
to be marginally consistent both with Planck Collaboration 
et al. (2015a) and the CFHTLenS analysis by Heymans et al. 
(2013). Because the constraints on as and flm are degener¬ 
ate, we combine both parameters into the projected param¬ 
eter Ss = (Ts(f2m/0.3)°'®. For the baseline ACDM model we 
obtain a best-fitting value of Ss = O.TOSIo'qsI- Marginal¬ 
ization over a photometric redshift bias per tomographic 
bin increases the errorbars on by ~ 7 per cent. Further¬ 
more, we compare our constraints on cosmological param¬ 
eters with other CFHTLenS studies and the recent result 
from DES and we find general agreement. Combining the 
shear likelihood with the Planck likelihood (TT+lowP) and 
sampling a simple six parameter ACDM model breaks the 
degeneracy between flm and erg and yields the constraints 
flm = 0.300 ± 0.011 and erg = 0.818 ± 0.013. These con¬ 
straints are consistent with the ones derived from resam¬ 
pling the Planck-only likelihood and the errorbars decrease 
by « 19 per cent for Dm and « 10 per cent for erg. 

Data from larger weak lensing surveys such as the Kilo- 
Degree Survey (KiDS^^, de Jong et al. 2013, 2015; Kuijken 
et al. 2015), the Subaru Hyper SuprimeCam lensing survey 
(HSC^'‘), and the Dark Energy Survey (DES^^ Flaugher 
2005; Jarvis et al. 2015; Becker et al. 2015) are building 
up right now, and these surveys will reach full coverage in 
the next years. This development will culminate in the sur¬ 
veys carried out by the Large Synoptic Survey Telescope 
(LSST^®, Ivezic et al. 2008) and the spaceborne Euclid}^- 
survey (Laureijs et al. 2011). Given these surveys, we con¬ 
sider our analysis also as a proof of concept in preparation 
for the order(s) of magnitude increase in survey area, which 
also implies a significant reduction in statistical uncertain¬ 
ties. 
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APPENDIX A: INDICES AND DERIVATIVES 

In Section 3 we described the generalization of the quadratic 
estimator to include tomographic bins. This requires a great 
amount of indices in a strict notation. For brevity we 
switched rather quickly to a set of super-indices and we also 
refrained from showing the explicit forms of certain matri¬ 
ces. Here, we give now the explicit forms of these matrices 
and we also calculate the derivatives used, for example, in 
equations 11 or 15 in index notation. 

First, we start with specifying the components of the 
vector B which contains all band powers /3 of type 6 for 
each unique redshift bin correlation ^ as Note that 

the total number of unique correlations between redshift 
bins is ricorr = n^iriz + l)/2, because all cross-correlations 
contain the same information by construction. Likewise we 
define the components of the tensor G which encodes all 
geometric information of the field depending on the band 
power /3, the band type 9 and the redshift bin correlation 
as 


G(ep(iJ,u)(ab)(ij) — 




d/3) 


Ai 


„i„(/3) 2(£-|-1) 


(Al) 


We also note here that each realization of G for a given band 
power P of type 9 and correlation C, can be represented as a 
matrix G^e/s- We can write out the matrices I® and Q® for 
the EE-, BB-, and EB-band powers as (Hu & White 2001): 


jEE 

jBB 

jEB 

and 


f Jo C 4 J 4 , 
\ S 4 J 4 


S 4 J 4 

Jo — C 4 J 4 


j Jo — C 4 J 4 
—S 4 J 4 


—S 4 J 4 
Jo C 4 J 4 


f — 2 S 4 J 4 2C4J4\ 

‘1C4J4 2S4J4J ’ 


Q™ 

Q™ 


Jo J- 2C4 J 4 + Cg Js SgJg 

SgJg — Jo J- 2 C 4 J 4 — Cg Jg 

— Jo + 2C4 J 4 — Cg Jg —SgJg 

— SgJg Jo J-2C4 J 4 + Cg Jg 

—2sgJs 2Jo + 2c 8 Jg\ 

2Jo + 2csJs, 2sg Jg J 


(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 


Section 4). Moreover, we employ the shorthand notations 
Cn = cos{ni^) and Sn = sin(n(p), where ip is the angle be¬ 
tween the x-axis and the distance vector Vij between pixels 
i, j (cf. Section 4). Note that in equation A6 we corrected 
the misprint in the original reference pointed out by Lin 
et al. (2012). Note also, that each block in the matrices of 
equations A2 to A7 defines again a matrix in the indices i, 

j- 

The matrices Mj in equation Al map between the red¬ 
shift bins and their unique correlations. In order to construct 
them, we start with the standard basis for fixu matrices 
with fijf & (l,...,nz). For example, the standard basis for 
riz = 2 can be written explicitly as: 

= (0 0) ’ ^ (0 0) ’ 

621 = , 622 = (^° J) . (A9) 

The index pairs (/i, v) can be trivially mapped to only one 
index which yields for the example above, i.e., p,,v £ 
( 1 , 2 ): 

(1,1) ^ 1, (1, 2) ^ 2, (2,1) ^ 3, (2, 2) ^ 4. (AlO) 


Imposing now, however, the symmetry condition (^, v) = 
which guarantees that for Uz redshift bins we only 
consider ricorr = nz(nz + l)/2 independent correlations, 
yields the symmetric mapping matrices: 


Ml 


M 2 


Mg 



= 612 J- 021 , 



(All) 

(A12) 

(A13) 


which implies the following mapping from (/r, zz) to 

(1,1) ^ 1, (1, 2) = (2,1) ^ 2, (2, 2) ^ 3. (A14) 

Next we construct the signal matrix C®'® as the sum over 
bands /?, band types 9, and redshift-correlations C, of the 
product of the band power vector with the geometry 

matrices Gpeq, 

E BcepG^gp(fzi.)(ab)(if) ■ (A15) 

C.0,/3 

Note that the full covariance matrix C also includes contri¬ 
butions from the shape noise matrix (cf. equation 9) 

which is constant. Thus if we wish to take the derivative 
of the full covariance matrix with respect to every possible 
band power combination B^fziz)(pe), this constant noise term 


vanishes and 

we are left with 


^G(izv)(ab)(ij) 


(A16) 

dBcep 

dBcep 


~ Gc9P{iJ,v){ab){ij) 

= Dce/3 = Da ■ 

(A17) 


In these equations we suppressed the argument of the Bessel 
functions J„ which in each case is the product ivij, where 
rij = |ni — rijl is the distance between pixels i, j (cf. 


In order to simplify our notation with respect to the Newton- 
Raphson algorithm we introduced in the last step the super¬ 
index A: each specific index combination {C,9P) can be 
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mapped to a single index i.e., we denote a specific 

derivative matrix now as Da instead of Hence the 

components of the generalized Fisher matrix F can be writ¬ 
ten as: 

Fab = I TriC-^DAC-^Dfl]. (A18) 

All other equations employed in the Newton-Raphson algo¬ 
rithm still hold with respect to this new set of super-indices 
iA,B). 

Finally, it only remains to write out explicitly the 
derivatives of the full covariance matrix C with respect to 
the power at an integer multipole 1. This is required for the 
calculation of the window function matrix (cf. equation 14) 
in which the derivatives enter in computing the trace 
matrix T (cf. equation 15): 


0 B(l) 

(A19) 

2 {ab)(ij)] 

(A20) 

where we have used that 



(A21) 

-f- 2 (^ab)(ij)] • 



APPENDIX B: ADDITIONAL FIGURES 

In the following figures we show the extracted E- and B- 
modes for each CFHTLenS patch individually. Note that 
these E-mode signals enter directly in the likelihood sam¬ 
pling whereas the combined signal presented in Fig. 4 serves 
just for illustrative purposes. 


For example, consider again two redshift bins, i.e., f £ (1, 2, 3), 
from which we wish to extract four band powers, i.e., /3 G 
(1, 2, 3, 4), of a single band type EE = 6 = 0. Then we can map 
each unique combination of (^6/3 to an integer A G (0,..., 12). 
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Figure Bl. Measured E-mode band powers in tomographic bins for the CFHTLenS patch Wl. From left to right we show the auto¬ 
correlation signal of the low redshift bin (blue), the cross-correlation signal between the low and the high redshift bin (orange), and the 
auto-correlation signal of the high redshift bin (red). The low redshift bin contains objects with redshifts in the range 0.5 < 21 ^ 0.85 
and the high redshift bin covers a range 0.85 < 22 ^ 1.3. The Icr-errors in the signal are derived from a run-to-run covariance over 184 
independent mock data fields (cf. Section 5.2) whereas the extend in ^-direction is the width of the band. Band powers in the shaded 
regions (grey) to the left and right of each panel are excluded from the cosmological analysis (cf. Fig. 5). The solid line (black) shows 
the power spectrum for the best fitting, five parameter ACDM model derived in the subsequent analysis (cf. Section 7 and Table 4). 
Note, however, that the band powers are centred at the naive £-bin centre and thus the convolution with the band window function is 
not taken into account in this plot, in contrast to the cosmological analysis. 



Figure B2. Same as Fig. Bl but for B-mode band powers. Note, however, the different scale (linear) and normalization used here with 
respect to Fig. Bl; for reference we also plot the best fitting E-mode power spectrum as solid line (black). We show the measured B-modes 
as (black) dots with Icr-errors derived from the inverse Fisher matrix. Based on these signals we define the shaded regions (grey) to the 
left and right of each panel (cf. Section 6 for details). E-mode band powers in these regions are excluded from the cosmological analysis 
(cf. Fig. 4). 
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Figure B3. Same as Fig. B1 but for CFHTLenS patch W2. Open symbols denote negative values plotted with Icr-errors centred on the 
absolute value. 



Figure B4. Same as Fig. B2 but for CFHTLenS patch W2. Open symbols denote negative values plotted with Icr-errors centred on the 
absolute value. 



Figure B5. Same as Fig. B1 but for CFHTLenS patch W3. Open symbols denote negative values plotted with Icr-errors centred on the 
absolute value. 
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Figure B6. Same as Fig. B2 but for CFHTLenS patch W3. Open symbols denote negative values plotted with Icr-errors centred on the 
absolute value. 



Figure B7. Same as Fig. B1 but for CFHTLenS patch W4. Open symbols denote negative values plotted with Icr-errors centred on the 
absolute value. 



Figure B8. Same as Fig. B2 but for CFHTLenS patch W4. Open symbols denote negative values plotted with Icr-errors centred on the 
absolute value. 
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